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A FREQUENCY INDEPENDENT TECHNIQUE FOR EXTRACTING THE 

RIGID -BODY MOTION FROM THE TOTAL MOTION 

OF A LARGE FLEXIBLE LAUNCH VEHICLE 

By James C. Howard 
Ames Research Center 


SUMMARY 


A method has Been devised for extracting rigid-body motion from combined 
rigid-body and flexible -body motion. The method involves the use of process- 
ing functions and requires that the number of sensing elements be equal to the 
number of modes of motion considered. Rejection of spurious flexible -body 
motion is accomplished by making the processing function associated with each 
sensing element a prescribed function of the modal slopes or the modal dis- 
placements, When the processing functions are allocated as prescribed, the 
information reaching an operator’s display panel, or the summing junction 
where error signals are generated, will not be contaminated by flexible - 
body motion. 

In the first part of the study, it is assumed that the modal data are 
accurately known. However, uncertainties in the distribution of mass, flex- 
ural stiffness, and shear stiffness may give rise to errors in the calculated 
modal data. Additional errors may result from inaccuracies in the mathemat- 
ical model used to compute the modal data. An investigation of the influence 
of modal errors on displays of an example launch vehicle led to the following 
conclusion. When the rigid-body pitch attitude is equal to, or greater than, 
the amplitude of the bending motion at the nose of the vehicle, the error in 
the pitch-attitude display never exceeds 1 6 percent provided the modal errors 
do not exceed either +50 percent or -50 percent. 

The effect of modal errors on the stability of the controlled system is 
such that certain combinations of errors tend to degrade the stability of at 
least one of the bending modes; whereas, other combinations of modal errors 
tend to enhance stability. The system used for stability analysis included 
the rigid-body mode and two bending modes. For modal errors of 10 percent or 
less, and those combinations of modal errors that tend to degrade stability, 
no instability occurred in either mode when the loop was closed with nominal 
gain. The results indicate, however, that the first bending mode is more 
sensitive to errors in the modal data than are the second and, presumably, 
higher bending modes. These conclusions are based on the assumption that the 
nominal modal data are used without modification. Examination of the data 
reveals that it is always possible to insure bending mode stability by a suit- 
able modification of the processing functions to account for an estimated 
magnitude of modal error. This is equivalent to the use of structural feed- 
back to stabilize the bending modes in the presence of errors in the modal 
data. 



IHTRODU CTION 


The task of manually controlling a large flexible booster is complicated 
by the manner in which structural feedback contaminates displays of vehicle 
motion • Likewise, the control of fully automatic launch vehicles is rendered 
more difficult by the presence of spurious bending signals at the summing 
junction where error signals are generated. To eliminate or mitigate flex- 
ibility effects, it is necessary to devise a method of extracting rigid-body 
motion or bending motion from the total motion of a flexible body. Convention- 
ally, this problem has been solved by the use of filters to exclude or attenu- 
ate all signals at bending mode frequencies. However, the advent of more flex- 
ible launch vehicle configurations has led to a reduction in the gap between 
the frequencies of the closed-loop rigid-body mode and the frequencies of the 
elastic modes. Reduction of this gap severely limits the use of conventional 
filters and requires the invention and development of new techniques for 
control. Research on this problem at Ames Research Center has been stimulated 
by the necessity of providing a human operator with displays of rigid-body 
motion which are not contaminated by spurious bending information. 

The primary objectives of the present study are: (l) to describe a fre- 

quency independent method for extracting the rigid-body motion from the total 
motion of a large flexible launch vehicle, ( 2 ) to determine the effects of 
errors in the modal data on displays of rigid-body attitude and attitude rate, 
and (3) to assess the effect of modal errors on closed-loop system stability. 

Other approaches to this problem are described in references 1, 2, and 3* 


SYMBOLS 


aj error coefficient 

rigid-body transfer function 
Di first bending mode transfer function 

D 2 second bending mode transfer function 

e error signal 

E Young's modulus of elasticity 

Fj generalized force 

F(s) filter transfer function 

G^ attitude sensor transfer function 

Gi(s) ith attitude sensor transfer function 
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rate sensor transfer function 
second, moment of area 

■fi 

generalized stiffness 
pitching -moment coefficient 
attitude gain 

attitude rate gain 

ratio of attitude rate gain to attitude gain 
generalized force coefficient for the first "bending mode 
engine inertia coefficient for the first lending mode 
generalized force coefficient for the second lending mode 
engine inertia coefficient for the second lending mode 
mass per unit length 
generalized mass 
lateral force distrilution 

function for processing the output of attitude and attitude- 
rate sensors 

function for processing accelerometer outputs in order to obtain 
rigid-lody pitching acceleration 

function for processing accelerometer outputs in order to obtain 
rigid-lody translational acceleration 

increment in Ppi(^) 

generalized coordinate 

increment in the pitch attitude 

pitch -attitude display error due to errors in jth modal slope 

column vector consisting of the rigid-lody pitch attitude and 
the n generalized coordinates 

complex frequency, a + joo 
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servo transfer function 
time 

distance measured along the longitudinal axis 
center -of -gravity location 

lending deflection at section x and at time t 

thrust vector deflection 

determinant 

rate sensor damping 

structural damping ratio 

bending angle at sensor location i 

bending angle at the nose 

commanded pitch attitude 

feedback signal 

output of an attitude sensor 

bending angle at nose in jth mode 

column vector of measured angular displacements 

real component of frequency vector 

value of the jth modal displacement function at location 
jth modal displacement function 
cofactor corresponding to 

matrix operator, the elements of which are the modal slopes at 
the sensor locations 

value of the jth modal slope function at location x^ 

jth modal slope function 

cofactor corresponding to 

imaginary component of frequency vector 

break frequency of first-order filter 



to j natural frequency of jth structural mode 

ixy^ rate sensor frequency 

C 0 g break frequency of first-order servo 

Subscripts 

A accelerated motion 

B magnitude of the subscripted variable resulting from bending 

eg center of gravity 

C commanded value of the subscripted variable 

i location at which subscripted function is evaluated 

3 denotes which bending mode the subscripted modal function describes 

n general term of a series of modal functions 

N magnitude of the subscripted variable at the nose 

P pitching motion 

R rate information 

T translation 

Superscript 

( ) derivative with respect to time 

METHOD 


The method proposed in the present report for extracting rigid-body motion 
from the total motion of a flexible body uses the series expansion of the 
motion of the elastic system in terms of the normal modes and the generalized 
coordinates. In the present context, a normal mode of free vibration is the 
spatial function or shape a vibrating beam assumes when it is oscillating at 
one of its natural frequencies. The generalized coordinates are the time 
functions or the functions that describe the variation of modal amplitude with 
time. The normal modes of free vibration of an elastic system can be computed 
when the distribution of mass, flexural stiffness, and shear stiffness are 
known. Hence, the modal functions can be obtained for any elastic system 
being considered. By equating measured sensor outputs to the motion described 
in terms of a series of normal modes and generalized coordinates, a system 
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of equations is obtained which can he solved for the rigid-body components of 
the motion, if the number of sensors equals the number of terms of the series. 
It is not anticipated, however, that more than a few terms of the series will 
he necessary, since high frequency bending signals can he eliminated, if 
necessary, by conventional filtering techniques. The solution of this system 
of equations is presented in the form of a series of products of processing 
functions and sensor outputs. The processing functions operate on the sensor 
outputs in such a way that the sum of the processed outputs gives the rigid- 
hody motion. 

The advantages of the proposed method are its simplicity and the absence 
of any need for on-hoard computational capability. All that is required is 
that the processing function associated with each sensor he a prescribed func- 
tion of the modal data at the sensor location. Since the modal data is a 
known function of position along the length of the vehicle structure, the 
processing functions may he precomputed, and hence no auxiliary on-hoard com- 
puting equipment is required. By the use of this method, separation and rejec- 
tion of bending motion is accomplished even in situations where a bending 
frequency coincides with the frequency of the closed-loop rigid-hody mode. It 
is not necessary to know the frequencies of the elastic modes since the method 
does not rely on the relative frequencies of the closed-loop rigid-hody control 
mode and the elastic modes. Instead, the method uses the known modal proper- 
ties of the structure being controlled, and consequently is well adapted to 
situations where a body bending frequency coincides with or differs only 
slightly from the control mode frequency. 


ANALYSIS 


Equations of Motion 

“When a beam is vibrating under the influence of a concentrated or a dis- 
tributed forcing function, the total displacement can be described in terms of 
the normal modes of free vibration and the generalized coordinates. A normal 
mode of vibration is the space function or shape the vibrating beam assumes 
when it is oscillating at one of its natural frequencies. In figure 1 are 
shown the first three bending modes of a uniform free -free beam. The general- 
ized coordinates are the functions which describe the variation of modal ampli- 
tude with time. These functions can be obtained by solving the partial 
differential equation of motion for an elastic beam which is moving under the 
influence of concentrated or distributed forcing functions. If shear deforma- 
tion and rotary inertia are neglected, the plane elastic motion of a continuous 
beam is described by the following partial differential equation, see 
reference 4. 
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3 2 y(x,t) 


a x 2 


+ m(x) 


5 2 y(x,t) 

at 2 


P(x,t) 


( 1 ) 
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The normal modes of vibration are obtained by solving the partial differential 
equation of motion for the free oscillations of an elastic beam. The relevant 
equation is 


bx. 2 


El(x) 


5 2 y(x,t) 

ax 2 


+ m 


(x> s , o 

5t 2 


( 2 ) 


A solution to this equation may be obtained by assuming that the variables are 
separable and by using the appropriate boundary conditions for the beam. A 
free -free beam is characterized by the absence of end constraints, that is, 
the bending moments and shear forces are zero at each end. With the use of 
these boundary conditions, equation ( 2 ) may be solved to obtain the spatial 
functions which may then be used to describe the displacement of a beam which 
is undergoing forced oscillations. In terms of the displacement and time 
functions, the total bending displacement is given by 


y(x,t) 


cpj(x)qj(t) 


(3) 


where cpj(x) are the modal displacement functions and qj(t) are the general- 
ized coordinates, or the functions which describe the variation of modal ampli- 
tude with time. 


Differential Equation for the Generalized Coordinates 

Upon substitution from equation ( 3 ) in equation (l), the following equa- 
tion is obtained 


5x 2 


El(x) 



X ^(x)qj(t) 
0=1 


Cj \ 1 

+ m(x) ^2 ) qVj(x)c[j(t) 


P(x,t) 


(M 


When each side of equation (4) is multiplied by cp^(x) and integrated over 
the length of the beam, taking advantage of the orthogonal property of the 
normal modes, the equation of motion for the jth mode assumes the following 
form 


Mjqij(t) + Kjqj(t) = Fj(t) 


(5) 


where Mj is the generalized mass of the beam in the jth mode of vibration. 
A dot denotes differentiation with respect to time. The generalized stiffness 
and the generalized forcing function in the jth mode of vibration are Kj 
and Fj, respectively. These quantities are defined as follows 
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( 6 ) 


M 


3 =jT m(x)cpj 2 (x)dx 


K J = . 


cpi(x) 


'4 El(x) tsM 

-C)x s OX^ . 


dx 


(7) 


Fj(t) = J' P(x,t)cpj(x)dx (8) 

where the integrals are taken over the length of the "beam. Equation ( 6 ) does 
not include the effects of rotary inertia, nor does equation (7) reflect the 
influence of shear deformation. Although shear deformation and rotary inertia 
should, in general, he included in any mathematical model which is being used 
to deter min e modal data, a discussion of these effects is not relevant to the 
present investigation, where the main objective is to process measured sensor 
outputs in such a way that rigid-body motion is extracted from combined rigid- 
body and bending motion. The preceding simplified presentation is used only 
to elucidate the arguments which follow. For a complete derivation of the 
equations of motion of an elastic beam, the reader is referred to references 
1 through 8. Equation ( 5 ) May "be rewritten as follows 


dj(t) + C0j 2 qj(t) 


Fj(t) 

M j 


(9) 


CO; 


El 

Mi 


( 10 ) 


where wj is the natural frequency of the jth free -free mode. The external 
forcing function consists of all aerodynamic forces, thrust forces, engine 
control-servo inertia forces, and propellant sloshing forces. 

Structural damping .- In practice, an elastic beam will possess some dis- 
sipative forces which provide damping. Since the dissipative energy is usually 
small in comparison to the elastic and kinetic energies, the lower modes will 
be very lightly damped. The dissipative force can be taken into account by 
a dding a small viscous damping term to equation (9)- When this is done equa- 
tion ( 9 ) assumes the following form 

q. + 2 CjC 0 jqj + coj 2 qj = ^ (ll) 

d 
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where £ j is the damping ratio in the jth bending mode, A value of 
- 0.005 "was assumed in the present study. 

Total displacement of an elastic beam which describes forced oscillations.- 
In terms of the normal modes and the corresponding generalized coordinates, 
the total displacement can be expressed as follow's 

n 

y(x,t) - cpj(x)<lj(t) (12) 

j=T,P,i 

where the normal modes are understood to include the rigid-body mode of trans- 
lation of the center of gravity, and the rigid-body mode of rotation of the 
beam about its center of gravity* In this formulation, the normal mode of 
translation is unity and the corresponding generalized coordinate is q T . In 
the rotational mode, the normal mode is x - x C g and the corresponding 
generalized coordinate is qp* Equation ( 9 ) can still be used to describe the 
rigid-body modes, provided 


— Wp — 0 


( 13 ) 


where cop is the frequency of the translational mode and oop is the frequency 
of the pitching mode. Note that all the aerodynamic forces are included in the 
external forcing function Fj(t). In view of these comments, equation (12) 
can be expanded as follows 


n 

y(x,t) = q/p + (x - x cg )q p cpj(x)q.j(t) (l4) 

j=i 


Separation of Rigid-Body and Bending Motion 

Angular displacements *- Let it be assumed that equation (2) has been 
solved and that modal displacement functions and modal slope functions are 
available. If 0 ^ denotes the measured output from an attitude sensor 
located at station i on the structure, then on differentiating equation (l4) 
with respect to x, the total angular displacement at location i is obtained 
in the following form 


9y dqpi(xi) 


+ 


dcpn( x i) 

<3x 


<Jn 


(15) 
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that is 


n 


e i 


= I'D + 


ScpjUi) 

c)x 


Clj(t) 


(16) 




If 6cpj(xi)/dx is denoted by tij* 
the jth bending mode at location i 
follows 


where ti < = ti(xi) is the modal slope for 

-LJ ’ J 

equation (16) may be rewritten as 


n 


e ± = <it 


^ij(x)qj(t) 


(IT) 




The measured output 0^ is seen to consist of (n+l) unknown components, that 
is, there is the unknown pitching angle q-p and the n generalized coordi- 
nates qj. However, a solution to equation (17) is possible if (n+l) sensors 
are used at (n+l) distinct locations. In this case there would be (n+l) 
measured outputs 0q, giving rise to (n+l) equations for the (n+l) unknowns as 
follows 
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• • • ^n+i,n/ 


(in / 


(18) 


This matrix equation can be solved to determine the unknown rigid-body rota- 
tion and the n generalized coordinates as functions of the measured outputs 
from the sensors and the known modal data. Equation (l8) may be written in 
abbreviated form as follows 


[©] = [\|r][Q] 


(19) 
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where [©] is a column vector of measured angular displacements, and [Q] is a 
column vector consisting of the rigid-body pitch attitude and the n general- 
ized coordinates. To simplify the formulation, the first column of the matrix 
operator [\|r] may he redefined to give 
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( 20 ) 


where \|^ 0 =• 1 for j-1, 2, . . . , n+1. From equation (l 9 ), the column vector 
[Q] is obtained in the following form 

[Q] - [^] _1 [©] (21) 

The inverse of the matrix [\| /] is given by 


Y10 

¥20 
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^n 2 

+n+i , 2 



• . . • « . 


\ ^2n ^3n * # # ^nn ^n+i,n j 

where Yj ± is the cofactor of the element ± in the matrix [\| /], and Ap 

-LJ - L «J 

is the determinant of [ij/] • Upon substitution from equation (22) in equation 
(21), the rigid-body rotation is obtained in the form of a series as follows 
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%> = (^io®i + ^2o®2 + • • • ^n+i,o®n+i^) 

that is, 

n+i 

q P = ^ X Yi ° 9i 

i=i 


( 23 ) 


The coefficient of 0q in equation (23) will he denoted by Ppi(t)> and- will 
subsequently be referred to as an attitude processing function, since this 
function is used to process the output from an attitude sensor located at xq. 
With this notation, equation (23) may be rewritten as follows 


where 


n+i 


%>= Y, P Pi (t)0 i 


1 = 1 


*io 
■pi' 1 T ' - An 


Eoi(t) = 


(2k) 


(25) 


Special case.- To illustrate the method, equation (18) will be solved for 
the case where i=3 and j=2. This is tantamount to the assumption that bend- 
ing modes higher than the second may be neglected. With these values of the 
subscripts, equation (l8) reduces to 
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In this case, equation (2^) consists of three terms. Hence, 
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where 
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Ap - [{^2X^32 +22+31 ) (Ul+32 " +12+31 ) + (+ 1 1 tpp ” t ip tpi ) ] 


(32) 


These are the attitude processing functions or gains to he used with the out- 
puts Q x , 0 2 , and 0 3 , respectively. When the outputs from the three attitude 
sensors are processed in the manner described, the information reaching an 
operator 1 s display panel, or the summing junction where error signals are 
generated will he rigid-body information only. Of course, this is contingent 
on the validity of the assumption that the motion of the system can he ade- 
quately represented by two bending modes. It can be seen that the attitude 
processing functions depend only on the modal slopes. Hence, these functions 
can be computed and stored prior to launch. Variation of modal data during 
flight may require that the appropriate processing functions be programmed to 
compensate for these variations. The sensitivity of the processing functions 
to changes in the modal data will be discussed later. 

Angular rates . - If 0 p denotes the output from an angular rate sensor 
located at x± on the elastic structure, then differentiation of equation (17) 
with respect to time yields the equation for the angular rates in the follow- 
ing form 

n 

o± = 4 P + t 1 j(x)4j(t) 

As in the case of angular displacements, the solution of the angular rate 
problem requires that (n+l) rate sensors be used at (n+l) distinct locations. 


(33) 
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The information obtained from this number of sensors is the amount required to 
permit separation of rigid-body rates from combined rigid-body and bending 
rates, by providing (n+l) equations for the (n+l) unknowns. The measured 
angular rates are related to the rigid-body pitching rate and the rates of 
change of the generalized coordinates in the same manner in which the angular 
displacements are related to the pitch attitude and the generalized coordi- 
nates. Hence, . . 

[®] - MCQ] (3*0 

where [©] is a column vector of measured angular rates, and [Q] is a column 
vector consisting of the rigid-body pitching rate and the rates of change of 
the generalized coordinates. The matrix [i|r] is defined in equation (20) where 
\|/q 0 = 1 for i=l, 2, . . . , n+l.. Equation (34) can be solved to obtain the 
column vector of unknown rates [Q] in terns of the measured rates and the 
known modal slopes 


[Q] = urt©] 


(35) 


where [\|r] is given by equation (22). It follows that the measured angular 
rates are processed in the same manner as the measured angular displacements. 
Therefore 


n+i 



i-i 


Ppi(+)<3i 


(36) 


where the processing functions Ppi(t) have the values given by equation (25)- 

The block diagram of figure 2 indicates how the rigid-body pitch rate may 
be extracted from the total motion of a flexible body. The three rate sensors 
are assumed to have transfer functions of the form shown. It is seen that the 
rigid-body pitch rate is obtained by weighting the output of each sensor by the 
appropriate processing function and summing. The functional forms Ppj_(t) of 

the processing functions are shown in the figure. The blocks Gi, i^l,2,3, 
represent the attitude sensors to be used in conjunction with the corresponding 
processing functions to give rigid -body pitch attitude. In cases where it is 
desirable to supplement the attitude and attitude -rate information with accel- 
eration feedback, the output of a suitable number of accelerometers must be 
processed in order to extract the rigid-body acceleration from the total 
acceleration of a flexible body. The method is similar to that described for 
obtaining rigid-body attitudes and attitude rates from the outputs of attitude 
and attitude-rate sensors. However, because of the fact that the rigid-body 
acceleration consists of a translational and rotational component, an addi- 
tional sensor is required to provide sufficient information for determining 
the rigid-body components of acceleration. The method is described in appendix 
A, where the processing functions required to extract rigid-body acceleration 
from accelerometer outputs are derived for the general case of the two rigid- 
body mode components and n bending mode components. 
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APPLICATION OF RESULTS OF ERROR ANALYSIS TO 
A LARGE FLEXIBLE BOOSTER 


Influence of Errors in the Modal Data 
on Displays of Vehicle Motion 

Since errors in the modal data give rise to errors in the processing 
functions, and hence to display errors and errors in the feedback control 
signal, it is necessary to examine the influence of these errors on displays 
of vehicle motion and on the stability of the closed-loop system. In this 
section, the influence of errors on displays of vehicle motion will be 
considered. 

The influence of errors in the modal data on displays of vehicle motion 
■was computed for a model vehicle. The model vehicle -which is representative 
of a typical large flexible space vehicle booster, and which was supplied by 
the Aero -Astro dynamics Laboratory of the Marshall Space Flight Center, is simi- 
lar to or larger than the Saturn vehicle, but is not intended to represent any 
specific vehicle. The modal data for this vehicle is shown in table 1 for a 
flight time of 80 seconds. To produce meaningful results, the locations chosen 
for the sensors represent possible sensor locations for the model booster con- 
figuration (see fig. 3)- For the following computations, the sensors were 
located at distances of 19-3* ^8.1, and 69*7 meters from the gimbal axis. 

Attitude errors and attitude -rate err ors resulting from inaccuracies in 
the modal data . - An IBM 7O9O digital computer was used to compute the attitude 
processing functions for a system with two bending modes. The initial computa- 
tions were for nominal values of the modal slopes at three distinct sensor 
locations. Since the influence of errors in the modal data on computed values 
of the processing functions depends on which modal slope or group of modal 
slopes is in error, the computations must be sufficiently comprehensive to 
include all possible combinations of modal slopes. For three attitude sensors 
and with the assumption that modal excitation is confined to the first two 
bending modes, a total of six modal slopes is required to compute the three 
processing functions. In this case, a total of (2 6 -l) combinations of modal 
slopes is possible. Hence, if the percent error assigned to each member of a 
group of modal slopes can assume a positive or a negative value, it would be 
necessary to compute each processing function 63 times for each percentage 
error assumed in the corresponding members of the group. However, if it is 
assumed that the percent error in each member of a group of modal slopes has 
the same sign, relationships between the errors in different groups of modal 
slopes can be established. Because of these relationships, the number of 
computations can be reduced as follows: If tij is assumed to have an error 
of n percent subject to the constraint that all other modal slopes are free 
from error, an error in the remaining slopes from the same mode can be related 
to the assumed error in tij* Ike slope at sensor location 1 in mode j is 
denoted by ^ x .. Subsequent to the introduction of an n percent error, the 
slope at x j becomes (l+n/lOO) If the jth modal slope function is 

now normalized to at sensor location 1, the slopes at sensor locations 

2 and 3, in mode j, assume the values \|r 2 j/(l+n/lOO) and \|r 3 j /(l+n/lOO) , 
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respectively. Hence, an error of n percent in is equivalent to an 

error of IT percent in the slopes j and i 3 j, where 


N = 100n 

100 + n 


( 37 ) 


Similar reasoning shows that an n percent error in i|r 2 j is equivalent to an 
N percent error in and t 3 j. Likewise, an error of n percent in 

is equivalent to an IF percent error in tij and Consequently, a given 

percentage error in a single modal slope has the same influence on the process- 
ing functions as an equivalent percentage error in each member of a group of 
two modal slopes as follows: 


til oc t2lt31 
tl2 a ^22^32 
t21 °c tllt31 
t22 or 1 12^32 
t31 ex fllt21 
t32 ex 1 12^22 


Therefore, out of a total of 15 pairs of modal slopes, the errors in 6 of 
these are related to the errors in the modal slopes taken one at a time. 

Hence, it is only necessary to consider errors in the remaining 9 pairs of 
modal slopes. As a consequence of similar reasoning applied to the remaining 
groups of modal slopes, it can he shown that out of a total of 63 possible 
combinations of modal slopes, it is only necessary to consider 24 of these. 

Of course, these computations would apply to only one set of sensor locations. 
If a range of sensor locations were considered, the computations would have to 
be repeated for each set of locations. 


Errors in the pro cessing functions interpreted as display _er_rors . - By 
referring to figure 4, it can be seen that the attitude error per degree of 
attitude sensor output, or the rate error per degree per second of rate sensor 
output, is simply the error A Pj_ induced in the processing function by 

errors in the modal data. The combined effect of these errors is given by 


3 

APpi^i - 
i=i i=a 

\[ f true value of \|/ 

A\[r assumed error in \| / 

APpi PpiU + Ai|r) - Ppi(i) 
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^ ^pi^p + ®Bi) 


(38) 



1 


is the bending angle at sensor location i, that is 


2 

°Bi * + *iA> * X ♦ijlj 


The equation for 5qp may be rewritten as follows 

3 3 

6<lp = QLp ^ ^pi + ^ %i ^pi 
i-i i-i 


o 

It can be shown that ^ APp^ = 0; therefore this equation reduces to the 


simpler form 


1-1 


3 3 2 

A AV . 

PI 


&q P “ ) ^Bi ^pi ~ } 2 j ^Pi 


l-i 


l-i j-i 


therefore 


&q P /_, 2, 4/1 j ^Pi “ ^ 

j=i 1=1 j=i 


(39) 


where 


a j 


\I . VP- 
y iJ 


Pi 


i=i 


(^0) 


This equation may be used to relate &qp to the bending angle at some refer- 
ence location. If the nose of the vehicle is chosen as the reference location, 
and if Q bn denotes the bending angle at this location, then 


^BN 


Vi^ 


8 hj 


J=1 


1 


where 0-^^ = ^jq.j is the bending angle at the nose in the jth mode. Equa- 
tion ( 39 ) may be rewritten in the following form 
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2 


2 


■where 


Bq.i 


•p 


C- 

-1 


+1 




NJ 






J=i 


♦.3 Hi 


5 %>j 0 Nj 


0=1 


Sq.i 


L P0 


_2_ 

%o 


(U) 


(^2) 


Numerical Results 

It is seen that &qpj is a measure of the error in the rigid-body pitch- 
attitude display per degree of bending at the nose in the jth mode. Like- 
wise, Stlpj may be interpreted as the error in the rigid-body pitch-rate dis- 
play, per degree per second of bending rate at the nose, in the jth mode. 

If there are no errors in the modal data, these coefficients are zero and the 
bending motion is completely filtered out. Wien the data corresponding to a 
flight time of 80 seconds was processed, the values of the coefficients shown 
in table 2 were obtained. The computed values of these coefficients are tab- 
ulated as a function of the percentage errors in a given modal slope or group 
of modal slopes. It is seen that when the amplitude of the rigid-body pitch 
attitude is greater than, or equal to, the bending motion at the nose of the 
vehicle, the error in the pitch attitude display never exceeds 1 6 percent if 
modal slope errors do not exceed either +50 percent or -50 percent. 


Influence of Errors in the Modal Data on the 
Stability of the Closed-Loop System 


Analysis . - From figure k it is seen that if 6 ^ denotes the total angu- 
lar displacement as seen at the ith sensor, then from equation (17) 


2 

g i “ 3-P + X * 13*3 

j=i 


(U) 


Therefore, 

~ Ip + + ^12 ^2 

^2 “ + ^21^1 + ^22^2 

£3 = Op + ^si^i + ^32 ^2 
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For the case considered, the air-frame dynamics may he represented by three 
parallel elements as shown in sketch (a). 



Sketch (a) 


Dr rigid-body transfer function 

Dx first bending mode transfer 
function 


D 2 second bending mode transfer 
function 


From equation (43) 

2 






m 


where & is the thrust vector deflec- 
tion. Therefore 


+ ^11^1 + tl2^2 
Os/b - fR + ^21^1 + ^ 22^2 

Q 3 /& = Dr + ^3lDl + ^32D2 


From figure 4, it is seen that if the rate feedback loop is omitted, the open- 
loop transfer function is given by 


^ = F(s)s(s)^G 1 (D R + | u Di + ti 2 D 2 ) 


P pi (\|r) + AP pl (\|/) 


+ G 2 (Dr + t 2 iDi + ^22D2) 


Ppi(^) t* APp^(\J/) 


+ ^"3 ( Dr + ^3lDl + t32D2) 


r pl 


(t) + Appiu) 


K r 


( 45 ) 


i = 1 , 2 , or 3 in first, second, or third group, respectively. 


feedback signal 
e error signal 

F(s) filter transfer function 

S(s) servo -transfer function 

Gp attitude sensor transfer function 

pitch attitude gain 
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With the assumption that the attitude sensors have the same transfer function, 
that is, G x = G 2 = G 3 = G A , equation (45) can he rewritten as 


-f - F(s)s(s)G A 


Br 


3 3 3 

^ ^pi(^) + -Pr ^ ^Fpi(^) + Di ^ ^iiFpx(^) 


L 1=1 


1=1 


1=1 


3 3 3 

+ Dl X APpi(t) + D 2 ^ ^i 2 Ppi(t) + Ds ^ ti 2 AP pi (t) 


1=1 


1=1 


1=1 


Kq r 


However, by using equations (29) through ( 32 ) and the definition of APpq(t) , 
it can be verified that 



i=i i=i i=i 


(46) 


and 


3 

y Ppi(t) - 1 ( 47 ) 

i=i 


(It is interesting to note that 


3 

^ AP pi (t) = 0. 
i=i 


This implies that in the 


absence of bending, the rigid-body display will not be affected by errors in 
the processing functions induced by errors in the modal data.) Therefore, 


d JL = F(s)s(s)G 


A 


D-^ + D] 




11 'S’piW 


+ 




i=i 


i=i 


K q.r 


As defined in equation (40) 

3 

ai = y \|ai AP pi (\|f) 
i=i 


3 [ 

a-2 = y ^ i2 AP pi (t) 

i=i 


(48) 


(49) 


(50) 
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(51) 


I 


therefore , 


^ = [F(s)S(s)G a (D r + a 1 D 1 + a 2 D 2 )]Kq p 


If attitude rate feedback is included, the open-loop transfer function must he 
modified accordingly. In this case 



e 


F(s)S(s) 




(G A + KrsGr) 



(52) 


where 


and 


K a 

K= - -2E 


R 


°R = 


^Rf 

s 2 + 2^R0JRS + COR 2 



+ 1 


Since principal structural and control system frequencies are very much less 
than ojj} it is permissible to assume for present purposes that Gr = 1. Like- 
wise, because of the absence of free gyro dynamics, the attitude transfer 
function also assumes the form G^ = 1. When these substitutions are made in 
equation ( 52 ), the open-loop transfer function assumes the form 




F(s)s(s)/lk + J ajDjW 


%s) 


K, 


dr 


J=1 


(53) 


Transfer function formulation ,- For the purpose of ascertaining the 
influence of the error coefficients aj on the stability of the closed -loop 
system, an approximate transfer function for the system shown in figure k was 
used. The filter transfer function F(s) and the servo -transfer function S(s) 
were both assumed to be first-order lags. In formulating the rigid-body trans- 
fer function Dr(s), engine inertia and aerodynamic effects were omitted. The 
first and second bending mode dynamics were treated as second-order systems 
in which engine inertia and structural damping were included. The open-loop 
transfer function given by equation (53) consists of the following elements: 
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F(s) = 
S(s) = 


oaf 

s + oaf 
k>S 

s + oa 


s 


Dn = ^ 
% S 2 


D, = 


Kn+ K 12 s 2 


1 s 2 + 2^ 1 oa 1 s + oa-j_ s 


K P H“ Kp P s 
s 2 + 2^ 2 oa 2 s + oa 2 2 


Dp - 


See table 3 for list of coefficients used. 

Error coefficients .- The error coefficients ai and a 2 are a measure of 
the influence of errors in the modal data on the dynamic response of the sys- 
tem. If there are no modal errors, these coefficients are zero and, conse- 
quently, the lending motion is completely filtered out. The determination of 
the extent of the influence of modal errors on the stability of the controlled 
system requires that equations (4-9) and (50) be evaluated. A digital computer 
program was formulated to evaluate these equations for a range of assumed 
errors in the modal slopes. The error coefficients which are tabulated as a 
function of the percentage errors in a pair of modal slopes in table 4 
are typical of the values obtained. The error coefficients corresponding to 
the errors in any other group of modal slopes are related by equation (42) to 
the Sq^ coefficients shown in table 2. In computing these coeff icients, it 
was assumed that the sensors were located at distances of 19- 3 j 48.1, and 69.7 
meters from the gimbal axis. A flight time corresponding to the condition of 
maximum dynamic pressure was also assumed. The influence of the error coef- 
ficients on the dynamic response of the system is given by equation (53) • 

This equation was evaluated by using the tabulated values of the error coef- 
ficients as input to a second digital computer program which was devised to 
provide frequency response and root-locus data. From the stability point of 
view, the effect of errors in the modal data may be thought of as changing a 
system with ideal (cancellation) compensation into one with unideal compensa- 
tion, the poles of the compensator being the bending mode poles. In the ideal 
case, when modal errors are absent, the poles and zeros of the compensator 
occupy the same locations in the complex plane. However, when modal errors 
are introduced, the zeros of the compensator move away from the poles. For 
example, if the modal errors are such that a given error coefficient aq 
assumes a negative value, the zeros of the corresponding compensator move away 
from the bending mode pole in the direction of increasing complex frequency. 
Alternatively, if the modal errors are such that an error coefficient assumes 
a positive value, the zeros of the corresponding compensator move closer to 
the origin of the complex plane. The separation of the poles and zeros of the 
compensator with modal errors introduces the possibility of an unstable 
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condition when the loop is closed. Examination of the tabulated error coeffi- 
cients reveals that for the sensor locations and flight time chosen, these 
coefficients can have the same or opposite sign, depending on which modal 
slope or group of modal slopes are in error; and also on whether the modal 
errors are positive or negative. Since the results contained in table k have 
all possible combinations of sign, these values of the error coefficients were 
used in conjunction with equation (53) to construct the root-locus plots shown 
in figures 5 through 18 . 

Error coefficients both positive . - In figures 5 through 9 it is seen that 
if there is a negative error in mode 1 at sensor location 1, and a negative 
error in mode 2 at sensor location 1, then both error coefficients are positive. 
In this case, both bending loci lie entirely in the left -half plane, indicating 
that both bending modes are stable for all magnitudes of modal errors. 

Error coefficients ax p ositive and a^ negative . - This combination of 
error coefficient signs is produced by a negative error in mode 1 at sensor 1, 
and a negative error in mode 2 at sensor 2. Figures 10 through 12 show that 
in these circumstances, the root -locus for the first bending mode lies entirely 
in the left-half plane, indicating that in the presence of such modal errors, 
the first bending mode is always stable. Although the locus for the second 
bending mode lies initially in the left -half plane, the figures show that as 
the error in the second mode is increased negatively at sensor 2, the stabil- 
ity of the second bending mode is degraded. Further increases of modal errors 
in the sense indicated may cause the second bending mode to become unstable. 

Error coefficients a x negative and a 2 positive .- A positive error in 
mode 1 at sensor 1 and a positive error in mode 2 at sensor 2 produces an error 
coefficient ax which is negative and an error coefficient a 2 which is posi- 
tive. As can be seen from figures I 3 through 15 the root-locus for the second 
bending mode lies entirely in the left -half plane, indicating that the second 
bending mode is always stable in the presence of such modal errors. In this 
case the locus for the first bending mode, which initially lies in the left- 
half plane, moves into the right-half plane as the modal errors are increased. 
Hence, modal errors of this kind can produce an instability in the first bend- 
ing mode. 

Error coefficients both negative . - If both the first and the second bend- 
ing modes have positive errors at sensor 1, both error coefficients assume 
negative values. Although both bending modes are stable for small errors in 
the modal data, the presence of modal errors of this type degrades the stabil- 
ity of both bending modes. If modal errors are sufficiently large, both bend- 
ing modes become unstable. Examination of figures 1 6 through 18, however, 
reveals that for the nominal range of gain considered the first bending mode 
is more sensitive to modal errors than the second. 

For the sensor locations and the flight time chosen, it is immediately 
evident that the system is stable for small percentage errors in the modal 
data. However, it is seen that if the error in a given modal slope or group 
of modal slopes is such that an error coefficient ai(i=l,2) assumes a posi- 
tive value, such modal errors have a stabilizing influence on the correspond- 
ing mode i; whereas, if the modal errors are such that an error coefficient 
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ai assumes a negative value, the errors will have a destabilizing influence 
on the corresponding mode i. This suggests that an attempt be made to bias 
the nominal values of the modal slopes in such a way as to insure that the 
error coefficients are always positive. Examination of the tabulated data 
reveals that it is possible to do this. For the case considered in this report, 
both error coefficients can be made positive by assigning a negative bias to 
modal slopes \|r i:L , \|f 2;L , an ^ a positive bias to the slopes \|/ 31 , ^ 22 ? axl ^ L 

^ 32 . This is equivalent to the use of structural feedback to stabilize the 
bending modes in the presence of errors in the modal data. 


CONCLUDING REMAEKS 


A method has been devised for extracting rigid-body motion from the total 
motion of a flexible body. The method, which does not rely on the relative 
frequencies of the closed-loop rigid-body control mode and the elastic modes 
as conventional filters do, is well adapted to situations where a flexible- 
body frequency coincides with, or differs only slightly, from the control mode 
frequency. 

An analysis of the effects of modal errors on motion displays of a typical 
launch vehicle indicates that when the amplitude of the rigid-body pitch atti- 
tude is greater than or equal to the amplitude of the bending motion at the 
nose of the vehicle, the error in the pitch-attitude display will always be 
less than 1 6 percent if the modal slope errors do not exceed either +50 percent 
or -50 percent. 

An analysis of the influence of modal errors on the stability of the 
closed-loop system indicates that certain combinations of modal errors tend to 
degrade stability of at least one of the bending modes; whereas, other combina- 
tions of modal errors tend to enhance bending mode stability. For modal errors 
of 10 percent or less, and those combinations of modal errors that tend to 
degrade stability, no instability occurred in either mode when the loop was 
closed with nominal gain. The results indicate, however, that for the nominal 
range of gain considered the first bending mode is more sensitive to errors in 
the modal data than the second and, presumably, higher bending modes. 

It has been shown that when the error coefficients are positive, stability 
is maintained without degradation. This suggests that an attempt be made to 
bias the nominal value of the modal slopes in such a way as to ensure that the 
error coefficients are always positive. Examination of the tabulated data 
reveals that it is possible to do this for the range of modal error considered 
in this report. This is equivalent to the use of structural feedback to stab- 
ilize the bending modes in the presence of errors in the modal data. 

It is not claimed that the results of the error analysis are complete. 

The assumption has been made that the errors in each group of modal slopes are 
always in the same direction, that is, all positive errors or all negative 
errors. It is more probable that one or more members of any group of modal 
slopes will be subject to errors in one direction? while the remaining slopes 
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of the group will have errors of opposite sign. However, the error analysis 
given is felt to provide a basis for the assessment of the effects of the 
modal error distribution. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., June 28, 1965 
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APPENDIX A 


ACCELERATIONS IN THE PITCH PLANE 


Let A± denote the output of an accelerometer which has its sensitive 
axis in the plane of the motion and perpendicular to the longitudinal axis of 
the vehicle. If second-order quantities are neglected, and if (xj_ - x cg ) is 
the distance of the ith accelerometer from the center of gravity of the 
vehicle, then the linear acceleration at the ith accelerometer may he 
obtained "by differentiating equation (lk) twice with respect to time 


(Al) 


A-i = 


q T + (x d 


x cg ) %> 




cp ij <1 «j 


where q-p is the component of rigid-body acceleration in the translational 
mode, and cp. . = cp.(xj) is the modal displacement in the jth mode at sensor 
location i. 


Because of the fact that the rigid-body acceleration has two components, 
the number of sensors required to provide sufficient information for determin- 
ing the rigid-body components and the bending components is (n+2). Hence, the 
equations to be solved are 


/ l 


11 (x x - xcg) 9 X1 9 12 ... 9 in \ 


IM 

1 K, 1 

• 


1 (x 2 - Xcg) 92 i 922 * • • 92n ' 

• • • » • 1 • • 


% 

• 

• 

- 

« • « • • • • • 
• • » » ... , 


• 

An 


1 (Xn - x cg) 9 nx 9ns • • • 9nn 


• 

An-f-i i 


i -L ( x n+i “ x cg) 9n+i,x 9n+i,2 • • * 9n+i,n , 


l-i 

, Aa+2 1 


\ 1 ( x n +2 - x cg) 9n+2,i 9n+2,2 • • • 9n+2,n j 

\ 


(A2) 


This matrix equation can he solved to determine the unknown rigid-hody accel- 
erations and the n bending accelerations as functions of the measured outputs 
from the sensors and known modal data. To render these equations more 
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manageable, it is convenient to rewrite them in abbreviated form as follows 


[A] = [cp][Q] 


(A3) 


where [A] is a column vector of measured accelerations, and [Q] is a column 
vector consisting of the translational acceleration, the pitching acceleration, 
and the n bending accelerations. To simplify the formulation, the first two 
columns of the matrix [cp] are redefined as follows 


where 


and 


1 'PiT 

9 x p 

*11 * * * 

'Pin 

cp 2T 

9 2 P 

cp 2 i . . • 

'Pan 

7 3 t 

'PaP 

cp 31 • ♦ • 

• • • • 

• ... 

'Pan 

• 

* 

• 

i 'Pn+ijT 

• 

'Pn+ijP 

• ... 

'Pn+iji * * * 

• 

'Pn+i, 

j'Pn+a^T 

'Pn+a.P 

'Pn+Sji * * ‘ 

c Pn+2 , 



<P iT = 1 


IP = (x i 

- x ) 

eg' 

i = 1^2, . . . 

,n+2 


(A4) 


The column vector of unknown acceleration components is obtained from equa- 
tion (A3) in the form 


[Q] = [cp] - 1 [A] 


(A5) 


where the inverse of the matrix [cp] 


is given hy 
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. -1 


a a 


$iT 

$2T 

• • • $n+i,T 

$n+ 2 ,T 

®iP 

<t 2P 

. . . $ n +i,P 

$n+ 2 ,P 

$11 

*21 

. . . $n+i,i 

$n+ 2 , 1 

$12 

$22 

. . . $n+i ,2 

$n+ 2,2 


■ in 


f 2n 


^n+i-n ®n+2,n 


(a6) 


where 4>ij is the cofactor of the corresponding element in the matrix [cp] and 
Aj^ is the determinant of [cp]. 

Rigid-body translational accelera tion. - Upon substitution from equa- 
tion (a 6) in equation (A5) , the rigid-body translational acceleration is 
obtained in the form of the series 




_ 1 _ 

A A 



iT A i + + 


that is 


ip 


n+2 



i=i 


®n+i,T^n+i 


+ $ 


n+2. 



(A7) 


The coefficient of Aq in equation (A7) will be denoted by Pqiq(cp) and is the 
function to be used in processing the output from an accelerometer located at 
xq on the flexible structure, in order to obtain the rigid-body translational 
acceleration. Equation (A7) shows that when the output from each accelerometer 
is processed in this way, the sum of the processed outputs gives the rigid- 
body translational acceleration (see fig. 19). In terms of the processing 
functions Pipq(cp) , equation (A7) may be rewritten as follows 


where 


n +2 


0-T 


= p Ti( c P)A i 

i=i 




£iT 

a A 


(AS) 


(A9) 
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Rigid - body pitching acceleration . - Equations (A5) and (A 6) give the rigid- 
body pitching acceleration in the form of the series 


Sp = El l $ iP A i + + 


+ ^n+i,P^4a^n+2,P^n-te 


that is 


It 


n+2 



1=1 


(A10) 


The coefficient of Aj_ in equation (A10) is the function to be used in proc- 
essing the output from an accelerometer located at in order to obtain the 

rigid-body pitching acceleration. This processing function will be denoted by 
Ppi(cp) • Equation (A10) shows that the sum of the processed outputs gives the 
rigid-body pitching acceleration. In terms of the processing functions 
P 1 • (cp), equation (A10) may be rewritten in the following form 

ir 

n+H 

= X p p jL ^ cp '* Ai 

i=i 

where 

p pi(<?) = ~ 


Similarly, by using equation (A 6) to obtain the appropriate processing func- 
tions^ it is possible to determine modal accelerations from which may be 
derived modal rates and modal displacements. These quantities may be used in 
a feedback loop to supplement attitude and attitude rate information. An 
interesting application which suggests itself is the use of structural feed- 
back to suppress or attenuate flexible body motion by augmenting the general- 
ized mass, the structural damping, and the structural stiffness. See 
figure 20. 
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TABDE I . - BENDING MODE DATA 




t - 80 sec 



meters 

*cp 1 (x) 

t x (x) 

*<P 2 (x) 

t 2 (x) 

0.1 

1.1219 

0.05065 

1.1585 

0.06504 

2.5 

1.0000 

.05097 

1.0000 

.06693 

4-9 

.8773 

.05122 

.8373 

.06837 

7-3 

.7543 

.05151 

.6728 

.06985 

9-7 

.6290 

.05283 

.4969 

.07623 

12.1' 

.5012 

.05346 

.3099 

.07843 

14.5 

.3738 

.05284 

.1254 

.07510 

16.9 

.2449 

.05386 

-.06336 

.07976 

19.3 

.1179 

.05192 

-.2431 

.06968 

21.7 

-.003921 

.04955 

-.3947 

.05402 

24.1 

-.1141 

.04493 

-.4790 

.02767 

26.5 

-.2197 

.04294 

-.5284 

.01334 

28.9 

-.3199 

.04058 

-.5429 

-.000884 

31-3 

-.4151 

.03867 

-.5324 

-.007722 

33-7 

-.5054 

.03662 

-.5063 

-.01393 

36.1 

-.5909 

.03604 

-.4647 

-.02850 

38.5 

-.6718 

.03127 

-.3838 

-.03867 

40.9 

-.7406 

.02625 

-.2836 

-.04463 

43*3 

-.7982 

.02171 

- . 1686 

-.05102 

45.7 

-.8446 

.01690 

-.03958 

-.05626 

48.1 

-.8763 

.008959 

.09996 

-.05935 

50.5 

-.8881 

.001829 

.2429 

-.06033 

52.9 

-.8864 

-.003258 

.3897 

-.06186 

55.5 

-.8727 

-.008119 

.5389 

-.06227 

57-7 

-.8476 

-.02174 

.6878 

-.06183 

60.1 

-.8090 

-.02364 

.8399 

-.07291 

62.5 

-.7397 

-.03391 

1.0081 

-.06786 

64.9 

-.6468 

-.04356 

1.1656 

-.06302 

67.3 

-.5299 

-.05408 

1.3085 

-.05554 

69.7 

-.3864 

-.06559 

1.5290 

-.04374 

72.1 

-.2158 

-.07636 

1.5155 

-.02809 

74.5 

-.02079 

-.08593 

1.5615 

-.009608 

76.9 

*1957 

-.09425 

1.5607 

.OIO32 

79-3 

.4306 

-.1013 

1.5119 

.O3O27 

81.7 

.6809 

-.1071 

1.4160 

.04946 

84.1 

I.OO37 

-.1430 

1.1845 

.1235 

86.5 

1.3616 

-.1551 

.8365 

.1638 

88.9 

1.7457 

-.1645 

.4044 

.1952 

91.3 

2.1338 

-.1595 

-.2044 

.2263 

93*7 

2.5204 

-.1624 

-.6694 

.2431 

96.1 

2.9174 

-.1697 

-1.23411 

.2730 

98.5 

3.3323 

-.1687 

- 1.9730 

.2763 

100.9 

3.7321 

-.1644 

-2.6196 

.2614 

103.3 

4.1197 

-.1487 

-3.2176 

.2373 


*Normalize& to unity at gimbal point . 



TABLE II . - DISPLAY ERRORS CALCULATED FOR GIVEN ERRORS IN MODAL SLOPES 


Percentage error 
in modal slopes 

Display error per degree 
of "bending at the nose of 
the vehicle, in mode j 

Percentage error 
in modal slopes 

Display error per degree 
of bending at the nose of 
the vehicle, in mode j 

*11 

5 Qpi 

j &qp2 

*12 

| s 9pi 

6q p2 

-50 

-O.O7568 



-50 


0 

0.08178 

-4o 

-.06005 



-40 



.06127 

-30 

-.041{67 



-30 



.04324 

-20 

-.02954- 



-20 



.02720 

-10 

-.01465 



-10 



.01288 

0 

0 



0 



0 

10 

.01442 



10 



-.01164 

20 

.02863 



20 



-.02221 

30 

.04261 



30 



-.O3I86 

4o 

.05638 



40 



-.04071 

50 

.0699^ 



50 



-.04880 

*21 

&q.pi 

| & q.p2 

*22 

i 5g Pi 

6qp2 

-50 

-.00708 



-50 



-.04227 

-40 

-.00561 



-40 



-.03067 

-30 

-.00416 



-30 



-.02104 

-20 

-.00274 



-20 



-.01293 

-10 

-.00136 



-10 



-.00600 

0 

0 



0 



0 

10 

.00133 



10 



.00523 

20 

.00264 



20 



.00984 

30 

.OO392 



30 



.01393 

40 

.00518 



40 



.01759 

50 

.00642 



50 

\ 


.02087 

*31 

6qpi 

& q P 2 

*32 

sqpi 

6q p2 

-50 

.13490 

c 


-50 

c 


-.03041 

-40 

.09489 



-40 



-.02474 

-30 

.06327 



-30 



-.01888 

-20 

.03802 



-20 



-.01281 

-10 

.01731 



-10 



-.00652 

0 

0 



0 



0 

10 

-.01468 



10 



.OO676 

20 

-.02728 



20 



.01378 

30 

-.03822 



30 



.01981 

bo 

-.04781 



40 



.02865 

50 

-.05628 

\ 


50 

\ 


.03653 

*11*12 

6q.pi 

&q. P 2 

" * 11*22 

6qpi 

6qp2 

-50 

-.10262 

.08630 

-50 

-.05707 

-.05993 

-4o 

-.07586 

.06380 

-40 

-.04963 

-.04185 

-30 

-.05287 

.04446 

-30 

-.03968 

-.02717 

-20 

-.03292 

.02769 

-20 

-.02765 

-.01555 

-10 

-.01544 

.01298 

-10 

-.01426 

-.00662 

0 

0 

0 


0 

0 


0 

10 

.01374 

-.01155 

10 

.01470 

.00467 

20 

.02604 

-.02190 

20 

.02954 

.00771 

30 

.03712 

-.03122 

30 

.04426 

.00941 

40 

.04715 

-.03965 

40 

.05871 

.01001 

50 

.05627 

-.04732 

50 

.07275 

.OO973 

* 11*32 

6q p i 

6q p2 

* 12*31 

6q.pi 

"qps 

-50 

-.06960 

-.01946 

-50 

. 16261 

.07458 

-40 

-.05568 

-.01761 

-40 

. 10968 

.05764 

-30 

-.04194 

-.01478 

-30 

.07048 

.04158 

-20 

-.02821 

-.01094 

-20 

.04082 

.02660 

-10 

-.01429 

-.00604 

-10 

.01789 

.01277 

0 

0 

0 


0 

0 


0 

10 

.01486 

.00728 

10 

-.01420 

-.01173 

20 

.O3O51 

.01596 

20 

-.02560 

-.02250 

30 

.04722 

.02622 

30 

-.03481 

-.03246 

40 

.06530 

.03833 

40 

-.04237 

-.04163 

50 

.08514 

.05263 

50 

-.04855 

-.05010 

* 21*22 

6q.pi 

&qp2 

* 21*32 

6q p i 

bqps 

-50 

-.OIIO3 

-.04568. 

-50 

-.00955 

-.02891 

-bO 

-.00780 

-.03245 

-40 

-.00720 

-.02381 

-30 

-.00525 

-.02187 

-30 

-.00511 

-.OI837 

-20 

-.OO316 

-.01323 

-20 

-.OO316 

-.01260 

-10 

-.00146 

-.00607 

-10 

-.00147 

-.00649 

0 

0 

0 


0 

0 


0 

10 

.00128 

.00518 

10 

.00121 

.OO678 

20 

.00235 

.00969 

20 

.00222 

.01395 

30 

.OO33O 

.01361 

30 

.00289 

.02137 

40 

.00410 

.01707 

40 

.OO335 

.02916 

50 

.00484 

.02010 

50 

.00350 

.03721 



TABLE II.- DISPLAY ERRORS CALCULATED FOR GIVEN ERRORS IN MODAL SLOPES 

Concluded 


Percentage error 
in modal slopes 

Display error 
of bending at 
the vehicle, 

per degree 
the nose of 
in mode j 

'^ 22^31 

tqpi 

6q. P 2 

-50 

O.I3I7O 

- 0.00051 

-40 

.09206 

-.00995 

-30 

.06126 

-.01159 

-20 

.03705 

-.09418 

-10 

.01708 

-.OO523 

0 

0 

0 

10 

-.01493 

.00582 

20 

-.02818 

.01188 

30 

-.01015 

.01804 

to 

-.05097 

.02415 

50 

-.06093 

.03021 

^ 31^32 

&<ipi 

& %,2 

-50 

.11681 

-.01880 

-40 

.08534 

-.03569 

-30 

.05898 

-.02165 

-20 

.03645 

-.01521 

-10 

.01695 

-.OO712 

0 

0 

0 

10 

-.01493 

.00621 

20 

-.02818 

.01176 

30 

-.01001 

.01669 

10 

-.05057 

.02115 

50 

-.06019 

.02516 

'^11^12^22 

‘q.pi 

6Qp2 

-50 

-.08931 

.OI366 

-10 

-.06812 

.O326I 

-30 

-.01369 

.02305 

-20 

-.03134 

.01158 

-10 

-.01506 

.00695 

0 

0 

0 

10 

.01106 

-.OO636 

20 

.02717 

-.01226 

30 

.O39I8 

-.01770 

10 

.05104 

-.02280 

50 

.06200 

-.02756 

'I'll*' 12*32 

6q.pi 


-50 

-.09321 

.05608 

-10 

-.06967 

.03877 

-30 

-.04936 

.02528 

-20 

-.O3I34 

.01171 

-10 

-.01506 

.00611 

0 

0 

0 

10 

.01112 

-.00185 

20 

.02761 

-.OO83O 

30 

.01075 

-.01070 

1*0 

.05360 

-.01201 

50 

.06637 

_ -.01235 

11 * 21 ’’ 32 


lp2 

-50 

-.08299 

-.01736 

-10 

-.06530 

-.01639 

-30 

-.OI835 

-.01116 

-20 

-.03191 

-.01070 

-10 

-.01594 

-.00598 

0 

0 

0 

10 

.01594 

.00528 

20 

.03201 

.01606 

30 

.01855 

.02638 

10 

.06557 

.03839 

50 

.O83I6 

.05230 

hi* 22^32 

sq.pi 


-50 

-.01963 

-.08314 

-10 

-.01159 

-.06245 

-30 

-.03658 

-.04395 

-20 

-.02616 

-. 02/52 

-10 

-.01385 

-.01298 

0 

0 

0 

10 

.01511 

.01155 

20 

.03120 

.02387 

30 

.01808 

.O3IIO 

10 

.06557 

.O3936 

50 

.083I6 

.OI678 


Percentage error 
in modal slopes 

Display error 
of bending at 
the vehicle, 

per legree 
the nose of 
in mode j 

t22’'32 

&9pi 

tqp 2 

-50 

0 

- 0.08125 

-10 



-.06098 

-30 



-.04307 

-20 



-.02713 

-10 



-.01285 

0 



0 

10 



.OII63 

20 



.02225 

30 



.03194 

10 



.01083 

50 



.01901 

^ 11^12^21 


tq. P 2 

-50 

- 0.10793 

.08656 

-10 

-.O8OI3 

.06397 

-30 

-.O56I9 

.01154 

-20 

-.03544 

.02773 

-10 

-.01675 

.01298 

0 

0 


0 

10 

.01511 

-.01155 

20 

.02884 

-.02187 

30 

.04136 

-.03114 

40 

.05286 

-.03949 

50 

.06348 

-.04711 

4^1 1^12^31 

tq.pi 

& q P s 

-50 

.00013 

.08175 

-40 

.00215 

.06119 

-30 

.00262 

.04315 

-20 

.00222 

.02718 

-10 

.00128 

.01285 

0 

0 


0 

10 

-.00141 

-.01163 

20 

-.00289 

-.02225 

30 

-.00444 

-.03194 

1+0 

-.00592 

-.01083 

5 ° _ 

-.00710 

-.01901 

^11^2 1^22 

tq. &q. px 

&q.p2 

-50 

-.07169 

-.06I52 

-10 

-.06005 

-.01125 

-30 

-.01610 

-.02823 

-20 

-.03147 

-.01593 

-10 

-.01587 

-.00670 

0 

0 


0 

10 

.01580 

.0016I 

20 

. 0313 1 * 

.00759 

30 

.0 1610 

.00919 

40 

.06093 

.00973 

5 ° 

.07435 

. 009 J tO 

* 11 * 22*31 

■ ipi 

tips 

-50 

.01769 

-.03675 

-IO 

.01130 

-.02815 

-30 

. 00686 

-.01997 

-20 

.00377 

-.01256 

-10 

.00161 

-.00594 

0 

0 


0 

10 

-.00114 

.00527 

20 

-.00202 

.00999 

30 

-.00262 

.01420 

10 

-.OO3O9 

.01799 

50 

-.OO343 

.02141 

^11^31^32 

& q.pi 

tqpa 

-50 

.01715 

-.03312 

-10 

.01110 

-.02617 

-30 

.00686 

-.01955 

-20 

.00383 

-.OI3O6 

-10 

.00161 

-.00657 

0 

0 


0 

10 

-.00114 

.00670 

20 

-.00188 

.01365 

30 

-.00237 

.02082 

10 

-.00260 

.02828 

50 

-.00261 

.03603 
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TABLE III . - SYSTEM PARAMETERS 


K12 = 0.04412 (dimensionless) 
Kpp = 0.042116 (dimensionless) 
Uf = 17*5 (l/sec) 
w s = 30 (l/sec) 

K pi = O.85 (l/ sec 2 ) 

Kn = 26.355 (l/ sec 2 ) 

K21 = 23.96 (l/sec 2 ) 

101 = 5-039 (l/sec) 

102 = 11.99 (l/sec) 

t, = 0.005 (dimensionless) 

Kq.p =• O.75OO (dimensionless) 
k4 p = 0.8025 (sec) 


3 ^ 



TABLE IV.- VALUES OF ERROR COEFFICIENTS USED TO CONSTRUCT ROOT -LOCUS PLOTS 


Percentage error 
in modal slopes 

Error coefficients 

Percentage error 
in modal slopes 

Error coefficients 

\|r \|r 

Y 1 1 Y 12 

a i 

®2 

^11^22 

a i 


-50 

0.01526 

0.0201(8 

-50 

0.00849 

- 0.01422 

-40 

.01128 

.01514 

-40 

.OO738 

-.00993 

-30 

.OO786 

.01055 

-30 

.00590 

-.006447 

-20 

.001(895 

.00657 

-20 

.00411 

-.OO369 

-10 

.00229 

.OO3O8 

-10 

.00212 

-.001571 

0 

0 

0 

0 

0 

0 

10 

-.002043 

-.002741 

10 

-.002186 

.001108 

20 

-.00387 

-.00519 

20 

-.00439 

.OOI83 

30 

-.00552 

-.00741 

30 

-.OO658 

.002233 

40 

-.007011 

-.009409 

to 

-.OO873 

.002376 

50 

-.00837 

-.01123 

50 

-.01082 

.00231 
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UJ 

00 


Vehicle dynamics 


e = q P +5> j (x)q j (t) 



Figure 2.- Block diagram of system for processing attitude and attitude -rate information. 

















Figure 3 • - Model vehicle configuration showing possible sensor locations. 
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Figure 4. - Block diagram of the closed-loop system including the processing functions 





















<x 


i ,, x l 1 ]__ x^ 1 1 1 

32 28 24 20 " ~I6 12 8 
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Figure 8.- Effect of modal errors on root locus, both error coefficients 
positive, ai = 0.01128, a 2 = 0.01514. 










x Open -loop poles 
O Open -loop zeros 
□ Closed -loop poles 


i 



1st bending mode 








, A^ 22 =0.4f 22 
, A* 22 =0.3* 22 
, Af 22 =0.2f 22 

, a* 22 =o.i * 22 


Figure 15 . - Root -locus variations produced by positive error in mode 1 at 
sensor 1, and by positive error in mode 2 at sensor 2. 
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13.2 - 


x Open -loop poles 
O Open-loop zeros 
□ Closed -loop poles 



cr 


Figure 1 6 .- Effect of modal errors on the root locus, both error coefficients 
negative, ai = -0.002043^ &£ - -0.002741. 





14.0 - 




A *l2 =ai *l2 
A*, 2 = 0 . 2* (2 
Af , 2 = 0 . 3^12 
A ^| 2 = 0 . 4^ |2 


Figure 18.- Root -locus variations produced by positive 
sensor 1, and by positive error in mode 2 at 


error in mode 1 at 
sensor 1. 
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Figure 19 • - Block diagram of system for 










NASA -Langley, 1965 A“2040 


Airframe dynami 



Figure 20.- Block diagram of system for 


A = q T + (X X CG )q P +£<£.q. 


s ► 


^Al 2 

q T +(X| X CG )q P +<^|(X|)q|(t) 

S 2 +2^ a , cu a , S+cj a 2 




W A2 2 

q T +(X 2 -X C G)qp+<i()|(X2)qi(t) 

S 2 + 2 C A 2 2 S + ^2 




w A3 2 

q T +(X3-XcG)qp+^i(X 3 )q,(t) 

S 2 + 2 £ A3 oj A 3 S+cu A3 2 



rocessing accelerations in order to obtain the rigid-body 
hing acceleration. 
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